Heavy Ion Initial Conditions and Correlations Between Higher Moments in the 

Spatial Anisotropy 
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Fluctuations in the initial conditions for relativistic heavy ion collisions are proving to be crucial 
to understanding final state flow and jet quenching observables. The initial geometry has been 
parametrized in terms of moments in the spatial anisotropy (i.e. 62,63,64,65...), and it has been 
stated in multiple published articles that the vector directions of odd moments are uncorrelated 
with the even moments and the reaction place angle. In this article, we demonstrate that this is 
incorrect and that a substantial non-zero correlation exists between the even and odd moments in 
peripheral Au+Au collisions. These correlations persist for all centralities, though at a very small 
level for the 0-55% most central collisions. 



One proposal for modeling the initial geometry fluc- 
tuations in relativistic heavy ion collisions is utilizing a 
Monte Carlo Glauber calculation Using this model 
and the initial transverse positions of the struck nucleons, 
referred to as participants, one can calculate thepartici- 
pant eccentricity £2 and higher moments . In Q , these 
moments arc defined by: 



\J (r 2 cos(n4> P art)) 2 + {r 2 sin(n(j) part )} 2 



(r 2 ) 



(1) 



where n is the nth moment of the spatial anisotropy cal- 
culated relative to the mean position. The major axis 
associated with the nth moment is defined by: 



atan2((r 2 sin(n4> part )) , (r 2 cos(3(f> part ))) 



(2) 



We show an example event display in Figure Q] that 
includes the positions of the participant nucleons, and a 
visualization of the 62, £3, £4, and £5 moments. We have 
drawn the vector direction of each along the long- axis 
of the associated moment. The nth moment has an n- 
multiplct of directions that are equally valid, separated 
by 2ir/n. 



I. PREVIOUS RESULTS 

In , the authors state that the minor axis of triangu- 
larity (i.e. £3) is found to be uncorrelated with the minor 
axis of eccentricity (i.e. £2) in Monte Carlo Glauber cal- 
culations. This conclusion is repeated in 0] where it is 
stated that the orientation of triangular overlap shape 
due to fluctuations is random relative to the event-plane 
direction as determined by the elliptic anisotropy. This 
conclusion is then generalized in jj] where these authors 
state that the angles for the odd and even harmonics are 
uncorrelated. 



There have been many initial studies for how this 
spatial anisotropy translates into the final momentum 
space distribution of particles (e.g. @, d-Q). The fact 
that the angular orientations are uncorrelated between 
odd and even moments has an important impact on the 
methodology for experiments to determine these momen- 
tum anisotropy moments v n and for two-particle corre- 
lation measurements, with relevance for jet quenching 
observables. 
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FIG. 1. (Color online) Monte Carlo Glauber event display for 
a sample Au+Au collision. The grey circles represent the po- 
sitions of all nucleons. The green (red) circles are participant 
nucleons from the left (right) nucleus. The vector directions 
for the n = 2 — 5 and the spatial anisotropy pattern they rep- 
resent are overlaid. These are centered on the mean position 
as indicated by a black square. 
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FIG. 2. (Color online) Monte Carlo Glauber Au+Au distribu- 
tion for relative angular difference (£2 — £3) [radians] for differ- 
ent selections in N par t- A flat distribution at one from angular 
separations 0.0-7r/6 indicates totally uncorrelated quantities. 



II. METHODOLOGY 

We set out to confirm the findings of the above papers 
using the Monte Carlo Glauber framework. We have 
utilized the standard PHOBOS Monte Carlo Glauber 
code [loj with Woods-Saxon parameters and settings 
(Rq = 6.38 fm, a = 0.535 fm, d min — fm). We show 
the angular distribution between £2 and £3 for a set of 
Au+Au number of participant (N part ) selections in Fig- 
ure [5J Because of the two (three) fold symmetry for the 
e 2 (£3) moments, if the two angles are uncorrelated, the 
distribution should be flat at one from 0.0 - tt/6 radians. 
A clear correlation is found between the two moments in 
peripheral Au+Au events. The correlation strength then 
decreases for more central events [ill ]. 

Previous studies have noted that correlations exist be- 
tween the positions of participating nucleons and are in- 
trinsic to all Monte Carlo Glauber calculations [l2j . And 
at the same time, previous publications assumed that 
this would not lead to correlations in the vector direc- 
tion of odd and event moments since the odd moments 
are purely the result of fluctuations. 

We quantify the degree of correlation by taking the 
root-mean-square (RMS) of the distribution. If there 
were no correlation between these moments, the RMS 
of the angular separation would be 0.302 radians. The 
results as a function of N part are shown in Figure [3] For 
our initial view, we have zoomed in for N part < 100. The 
results confirm the strong correlation between the two 
angles (i.e. the large downward deviation in the RMS 
from the flat case) for N part ~ 10, the deviation then 
weakens for larger N part values. In this figure, the cor- 
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FIG. 3. (Color online) The root-mean-square (RMS) of the 
angular difference between £2 and £3 as a function of num- 
ber of participant nucleons, N par t, for Au+Au, Cu+Cu, and 
Si+Si collisions, as well as the expectation of no correlation. 
We plot Npart < 100 to highlight the significant correlations 
that are present at small N part . 
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FIG. 4. (Color online) The root-mean-square (RMS) of the 
angular difference between £2 and £3 as a function of num- 
ber of participant nucleons, N par t, for Au+Au, Cu+Cu, and 
Si+Si collisions, as well as the expectation of no correlation. 
We zoom in on the vertical axis to highlight the small devia- 
tions from no correlation that remain at large N par t. 



relation appears to disappear for N part > 50. However, 
in Figure 2J we show the full N part range zooming the 
vertical axis around the default value. We observe a very 
small remaining anti-correlation for N part ~ 80 — 150, 
which translates into an approximate 1% lower proba- 
bility of having the two moments within 0.0-7r/12 and a 
correspondingly higher probability of having the two mo- 
ments separated by 7r/12-7r/6. There is a similar magni- 
tude positive correlation for N part > 200. 

For the lower N part region where the alignment is 
strongest, the correlation may be due to small number 
fluctuations of the particular geometry in that impact pa- 
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Monte Carlo Glauber Results 
Au+Au @ 200 GeV 
No correlation case. 
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FIG. 5. (Color online) The RMS angular difference between 
the reaction plane, ^rp, and £3 as a function of N par t for 
Au+Au collisions. A weak correlation is found between the 
3rd-order moment and the reaction plane direction for small 
Npart- Note the zoomed vertical scale. 




Monte Carlo Glauber Results 
Au+Au @ 200 GeV 
Cu+Cu @ 200 GeV 
Si+Si @ 200 GeV 
No correlation case. 
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FIG. 7. (Color online) The root-mean-square (RMS) of the 
angular difference between £2 and £4 as a function of num- 
ber of participant nucleons, N par t, for Au+Au, Cu+Cu, and 
Si+Si collisions, as well as the expectation of no correlation. 
The result shows the same procedure applied to the differ- 
ence between the lowest even-order moments where a more 
significant correlation results from the average geometry. 
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FIG. 6. (Color online) The root-mean-square (RMS) of the 
angular difference between £2 and £5 as a function of num- 
ber of participant nucleons, N par t, for Au+Au, Cu+Cu, and 
Si+Si collisions, as well as the expectation of no correlation. 
The result demonstrates that a similar behavior holds for cor- 
relation between £2 and higher order odd moments. 



rameter range for Au+Au. Thus, also shown in Figure [3] 
are calculations from Si+Si and Cu+Cu collisions. One 
can see that the angular correlation largely tracks with 
Np ar t and thus it is the number fluctuations and not the 
average geometry that dominate the correlation. How- 
ever, the Si, Cu, Au results do not scale perfectly versus 
Npart j and thus the particular geometric configuration 
plays some role, as one might expect since the geometry 
correlates with the magnitude of these moments. 

Furthermore, the authors of [HI state that the fluctua- 
tions are random with respect to the reaction plane (the 
plane defined the line between the centers of the nuclei 
and the longitudinal axis). We have tested this as well 



and show the resulting RMS angular separation between 
the 63 and the reaction plane in Figure [5] We find a sim- 
ilar trend as shown above. The correlation between £3 
and the reaction plane is also restricted to small num- 
bers of participating nucleons. The correlation here is 
much weaker than what was found between £2 and £3. 

We also show the angular correlation between £2 and 
£5 in Figure [6l One again sees a strong correlation in pe- 
ripheral collisions and then smaller correlations for larger 
numbers of participants. Finally, in Figure[7l we show the 
angular correlation between £2 and £4. These two even 
moments are expected to be highly correlated because 
the initial overlap in mid-central collisions has an ap- 
proximately elliptical shape, which can be well described 
by a combination of aligned £2 and £4 moments if it has 
a large enough eccentricity. This correlation tracks with 
the geometry (i.e. impact parameter), and so the result- 
ing RMS values for Au+Au, Cu+Cu, and Si+Si do not 
track each other when plotted as a function of N part . 

In the case of nearly ideal hydrodynamics, the spatial 
eccentricities may be translated into directly measurable 
momentum anisotropies (v n ). The relative degree of cor- 
relation between vi and W4 is of great interest because it 
directly impacts a possible difference in the observed 
determined via the event plane method (using the sec- 
ond moment to determine the plane) and a two-particle 
Fourier decomposition method [lij]. The V3 moment is 
identically zero by symmetry when measured with re- 
spect to the second moment event plane. However, the 
vector direction correlation can be determined via event- 
by-event separate second and third-moment event plane 
measurements. Such measurements would provide po- 
tentially sensitive tests of the fluctuations in the Monte 
Carlo Glauber geometry. Despite the effect being largest 
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in peripheral events where nearly inviscid hydrodynamics 
may not apply, it should prove fruitful to see if such angle 
correlations persist (which do not require a linear trans- 
lation of e„ moments into v n momentum anisotropies) . 

III. SUMMARY 

In summary, we find that contrary to previous publica- 
tions, there is a definite non-zero correlation within the 
Monte Carlo Glauber calculation between the angular di- 
rections of £2 and €3 (and more generally between even 
and odd moments). A much weaker correlation between 
the reaction plane and €3 is also found. The effects are 
strongest for peripheral Au+Au events, and the measure- 
ment of such correlations may provide an interesting test 
of geometric fluctuations within the Monte Carlo Glauber 



calculation. It will be instructive to measure these possi- 
ble correlations and care must be taken to quote the exact 
sensitivity of such measures and extend them to the very 
most peripheral events. The multi-dimensional correla- 
tion between the magnitude of the eccentricity orders and 
their angular orientations may also prove important and 
the direct integration of Monte Carlo calculations of ini- 
tial state geometries may be warranted in many studies 
to account for the full set of correlations. 
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